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ABSTRACT 

A minimal particle content supersymmetric model with a discrete Z% symmetry, allowing lepton 
number violating terms, is studied. Within this model, the neutrino masses and mixing can be 
generated by lepton number violating couplings. Choosing parameters which correctly describe both 
the masses and mixing in the neutrino sector, we consider their repercussions in flavour violating 
radiative lepton decays, I — > l lr y. Such decays have not been observed and, accordingly, soon to be 
improved, upper bounds on their branching ratios exist. We do not assume dominance of either the 
bilinear or trilinear couplings. We note that certain parameter sets, which correctly describe the 
neutrino sector, will also generate observable branching ratios, some of which are already precluded, 
and suggest four such sets as Benchmark scenarios. 
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1 Introduction 



Processes which do not conserve lepton flavour, the flavour oscillations in the neutrino sector, 
have been observed [1]. This is in contrast with the charged sector, where no such observation 
has been made. The decays r — > ^7, cy and \i — > cy can be driven solely by the known lepton 
flavour violation in the neutral sector, the branching ratio will be small [2], however, well below 
current experimental limits, due to the magnitude of the neutrino mass. Noticing that in many 
extensions of the Standard Model this branching ratio increases greatly, 1 together with the 
fact that experimental bounds for these decays will soon be improved by several orders of 
magnitude, suggests that these decays are a valuable place to scrutinise the Standard Model 
and test theories which extend it. 

The model considered here, the Lepton Number Violating Minimal Supersymmetric Stan- 
dard Model (JZ-MSSM), is a well motivated extension of the Standard Model. When con- 
structing a supersymmetric model dangerous operators arise which, unless highly suppressed, 
would give rise to an unacceptable rate for proton decay. As such, a discrete symmetry must 
be imposed when constructing the Lagrangian [10,11]. We choose a discrete symmetry which 
allows lepton number violating terms, but not baryon number violating terms, resulting in 
proton stability. The lepton number violating terms violate R-parity. The superpotential for 
the t'-MSSM is given by 

Wl/_MSSM = tab {^^aPk^a^pEk + ^'ajk^aQ^ Afc,x + (Yj) jkQ'j'* '' H%Uk,x ~ ^a^H^j , 

(1.1) 

where Q°; x , Df, Uf , E i: H±, are the chiral superfield particle content, i = 1,2,3 is 

a generation index, x = 1,2,3 and a = 1,2 are SU(3) and SU(2) gauge indices, respectively. 

££= 3 = (iff, £"=1 2,3) an d hence, Aoj/c are the lepton Yukawa couplings and are lepton 

1 The effect of lepton flavour non-conservation from the charged slepton mass matrix in supersymmetric 
extensions of the Standard Model where a seesaw mechanism results in light Majorana neutrinos is noted in 
Ref. [3], this work is extended in Ref [4], where bounds for off-diagonal terms are calculated. In Ref. [5] the 
results for this model are correlated with neutrino masses and the (g^ — 2) data and in Ref. [6] the possibility 
of discriminating between different supersymmetric seesaw models is investigated. A bottom-up approach is 
considered in Ref. [7], resulting in predictions for the \i — > branching ratio. Methods for discerning models 
with heavy right handed neutrinos from R-parity violating models using a number of decays are studied in 
Ref. [8]. Renormalisation group effects due to R-parity violating couplings and their effect on the /1 — > 
branching ratio are considered in Ref. [9]. 
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number violating parameters, whose role will be considered later. [x a is the generalised dimen- 
sionful /x-parameter, with \x§ and Hi, i = 1, ...3 the lepton number conserving and violating 
parts respectively. Similarly, X' j k are the Yukawa couplings for the down-type quarks and 
A^ fc violate lepton number. Finally, (Yu)ij are Yukawa matrices with e a b being the totally 
anti-symmetric tensor ei 2 = +1. The are further sources of lepton number violation in the 
soft breaking sector, however, we present here merely the terms which will play a role in this 
paper, 

£ssb ^ B a u La h° 2 - B a ~e La ht + H. c. , (1.2) 

where B a is the four-component bilinear term B a = (_B ,-Bj) are Bi are the lepton number 
violating components. The full supersymmetry breaking part of the Lagrangian and the mass 
matrices along with further details of the model are presented in Refs. [12,13]. The calculation 
is performed in a basis where the vacuum expectation values for the sneutrinos are zero; a 
method for moving to this basis from the most general scalar potential and properties of this 
basis are outlined in Ref. [14]. 

A particularly noteworthy feature is that the current experimental values of neutrino mass 
squared differences and mixing can be accommodated within the model, being determined by 
the value of lepton number violating couplings in either the superpotential, or the supersym- 
metry breaking terms of the Lagrangian [13]. In fact it has been shown that, one, and only 
one, neutrino mass can be generated at tree-level [15-18] with the masses of the two remaining 
generations arising through radiative corrections, producing the hierarchy between solar and 
atmospheric mass differences in a convenient fashion. Furthermore, this scenario can arise 
from a Froggatt-Nielsen model in which the discrete symmetry is due to the breaking of a 
U(l) x [19]. 

Crucially, the operators which give rise to neutrino masses in this model may also give rise 
to lepton flavour violation in the charged sector. In this paper, we shall consider combinations 
of lepton number violating parameters that correctly reproduce the observations made in 
oscillation experiments. For these sets of parameters we shall investigate whether they would 
result in branching ratios of rare leptonic decays which would already have been observed, 
or would be observed by forthcoming experimental studies. If the rare leptonic decays are 
not observed, the improving bound will be valuable in precluding certain scenarios. We select 
scenarios in which the off-diagonal terms in the supersymmetry breaking scalar mass matrices 
are zero. It is possible, of course, even in the R-parity conserving MSSM, that this is not 
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the case and that these terms will lead to large branching ratios for lepton flavour violating 
decays [4]. The aim of this study is to examine, specifically, the effects of lepton number 
violating terms in the Lagrangian and the interplay between the charged and neutral sector. 
As such, we will examine the scenarios which are only present in the L'-MSSM and will not 
examine phenomena which have their origin in the R-parity conserving part of the Lagrangian. 

2 Experimental Results, Bounds and Prospects 

The results from oscillation experiments combine to describe the mass squared differences 
and mixing angles of the neutrino sector increasing accurately. The current 3a allowed ranges 
are [1] 

sin 2 12 = 0.24 - 0.40 , sin 2 6 23 = 0.34 - 0.68 , sin 2 13 ^ 0.041 , (2.3) 
Am 2 x = (7.1 - 8.9) x 10" 5 eV 2 , |Am^| = (1.9 - 3.2) x 10" 3 eV 2 . (2.4) 

In our analysis we choose Lagrangian parameters such that the neutrino mixing angles match 
the tri-bimaximal mixing scenario of Ref. [20], 

sin 2 12 = i, sin 2 23 = i, sin 2 13 = O. (2.5) 

The following bounds have been set on the branching ratios of fi — > e^y [21], r — > /ry [22] and 
t -» e7 [23]. 

Br(/i erf) < 1.2 x 10~ n at 90% CL (2.6) 
Br(r -» pry) < 6.8 x 10" 8 at 90% CL (2.7) 
Br(r e 7 ) < 1.1 x 10" 7 at 90% CL (2.8) 

Future experiments will probe these decays further. It is suggested [24-26] that the sen- 
sitivity to \i — > e7 will be improved to ~ 10 -13 ' -14 , and that the sensitivities to r — > ^7 and 
r — >■ e7 will reach ~ 10 _8,_9 . To conclude, we present the current experimental values for the 
branching ratios of r — > /i^ r ^, t — > ev T v e and fj, — > ez/ M z/ e [27], 

Br(r \lv t v^ = 0.1736 + 0.0005 (2.9) 
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Br(r — > ev T v e ) 
Br(/j -> evyp,,) 



0.1784 + 0.0005 
1 



(2.10) 
(2.11) 



3 Generic Diagrams for l % — ► L7 



At the level of one loop, three basic types of diagram contribute to the decay U — > 1^ and in 
each case, there is a fermion - boson loop. The external photon can be attached either to the 
fermion in the loop, the boson in the loop or the external leg (Fig. 1). The calculation, as 
shown in Fig. 1, was performed in Weyl notation. In this notation, the four-component spinor 



spinor, / 



where ei and e# are two-component left-handed spinors and the four-component 



fi 



denotes a generic fermion. The factors associated with the vertices are 



denoted by either Aik s or Bik s as shown in the diagrams, the charges of the particles are 
given by Q^ iip , the masses of the particles in the loop are m^ tV and the masses of the charged 
leptons on the external legs are given by m^j. For more information concerning calculations 
using two-component spinors, see Ref. [28]. Taking all possible combinations of arrows and 
neglecting diagrams with gauge bosons 2 it can be seen that, in agreement with Ref. [29,30], 
at leading order the branching ratio is given by 



48tt 2 



[|A L | 2 + |A fl | 2 ]r(Z i IjUii/j) , 



(3.12) 
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2 It can be seen that diagrams will be suppressed either by the magnitude of the neutrino mass or, in 
diagrams which contain lepton number violating operators, by the amount of mixing between the neutri- 
nos/neutralinos and charged leptons/charginos. 
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Figure 1: TTie ^nol state photon can be attached to either of the external fermions, the fermion 
in the loop, or the scalar in the loop. All combinations of helicities of fermions should be 
included in the calculation. 
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The contribution from each Feynman diagram was expanded under the assumption 
to 2 , to 2 « to 2 , to 2 . On-shell conditions could then be applied and terms proportional to 
TOj were neglected. The resulting expression can then be re-arranged to be seen to contribute 
to effective operators of the form 3 



£ eff 3 2A Rij e Rj a^e Li F^ + 2A Lij e Lj a^e Ri F liU . (3.13) 

Individually, diagrams produce terms contributing to different effective operators, but these 
cancel when all possible diagrams are considered. Following Ref. [29] and inserting the ex- 
perimental values given in Eqs. (2.9,2.10,2.11), it is then possible to move from this effective 
operator to a branching ratio for the rare decay, given by Eq. (3.12). 
3 We define o^ v = \ {a^a v - a v a^) and o» v = \(g^u v - d v a^) 
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Figure 2: The various possible combinations of particles which can be produced in the loop: 
chargino and neutral scalar; neutralino and charged scalar; down quark and up-type squark; 
up quark and down-type squark. 



4 Specific Diagrams for U — ► L7 



3 

The following combinations of particles can be produced inside the loop: chargino and neutral 
scalar; neutralino and charged scalar; quark and squark, as shown in Fig. 2. In the JZ-MSSM, 
mixing occurs between charged leptons/charginos and between neutrinos/neutralinos. For 
example, there are five charged fermions, 

*f-(3).«-i....... 

where i = 3, 4, 5 are the charged leptons e, \i and r. Similarly, there are seven neutral fermions, 

where i = 5, 6, 7 are the neutrinos. The two-component spinors comprising the quarks are 
denoted, 

di = | dLl \ , i = 1,...,3 , 
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Figure 3: Diagram contributing to I — > Z'7 £/ie only source of lepton number violating 
being the bilinear lepton number violating couplings in the superpotential, Hi. 




Each of the diagrams in Fig. 2 are in the same form as outlined in Sec. 3. The generic 
vertices and Bj^ s can be replaced by the appropriate Feynman rule, which are presented in 
the Appendix. The calculation is performed in the mass eigenbasis. The full mass matrices are 
diagonalised and the appropriate rotation matrices are calculated numerically and without 
approximation. In understanding the important physical contributions it is more useful, 
however, to present diagrams in the mass insertion approximation containing interaction state 
particles. The plots are based on a Fortran code which computes the full result. 

We will consider the role played by combinations of lepton number violating parameters 
by, first, investigating the case in which the bilinear lepton number violating coupling in the 
superpotential correctly produces the atmospheric mass difference (and the ratios between 
the three components ensure the mixing angles are reproduced correctly) and another, single, 
lepton number violating coupling sets the solar mass difference. Both sources of lepton num- 
ber violation will then combine to produce a diagram which contributes to a lepton flavour 
violating decay. Second, we will consider the scenario in which both the scale of the atmo- 
spheric mass difference and the solar scale are set by radiative corrections, and the bilinear 
lepton number violating parameters are set to zero. 

4.1 Atmospheric scale set by /Ai,2,3 

With only A*i,2,3 and all other lepton number violating couplings set to zero the atmo- 
spheric mass squared difference can be correctly reproduced; the solar mass squared difference 
is not generated and no observable branching ratios for / — > £'7 are generated. The non-zero 
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Figure 4: r — > /ji'y and Am 2 tm against fii - the dotted line (and right hand axis) show the 
magnitude of Amf tm with the light grey band indicating the current 3a allowed range; the full 
line (and left hand axis) show the branching ratio of ' r — > /try. 

^1,2,3 do bring about a branching ratio for I — > Z'7, through the diagram shown in Fig. 3. 
The fermion inside the loop is a mixture of the heavy neutralinos and the interaction state 
neutrinos. The amount of mixing between these interaction states is dependent on ^1,2,3, 
and also determines the mass of the tree level neutrino. The amount of mixing between the 
external leg interaction state charged higgsino and left handed charged leptons, is also deter- 
mined by /ii,2,3- As such, this diagram contributes to the U — > ^7 decay with branching ratio 
approximately given by, 



3 |A 0ii | 2 e 2 2 




T(k - IjViVj) . (4.14) 



We will consider this scenario. The /ii,2,3 parameter takes the values 

^ = ^| = ^= = 1.12 - 1.82MeV , (4.15) 

where ratios between components of fii are chosen such that the mixing angles in the PMNS 
matrix are generated to take the tri-bimaximal form, which are in agreement with current 
bounds, and are entirely generated in the neutral sector. A more detailed presentation of 
all the following scenarios is presented in Ref. [13]. The overall scale is varied and the re- 
sulting mass squared difference and branching ratio for r — > /ij are calculated and shown 
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Figure 5: Feynman diagrams generated when, firstly, fa induces mixing of interaction state 
neutrinos with gauginos and higgsinos and secondly, \jk ¥= generating neutrino masses 
radiatively. 



in Fig. 4. The light grey band (and right hand axis) shows the current 3a allow region for 
the atmospheric mass squared difference as given by Eq. (2.4). All other lepton number 
violating couplings are set to zero, and R-parity conserving parameters are fixed to be the 
SPSla benchmark point [31]. We note that, in agreement with Ref. [32], the branching ratios 
for li — > Zj7 are well below current experimental limits, as given in Eq. (2.7), and show the 
resulting branching ratio for r —> in Fig. 4. 



4.2 Atmospheric scale set by Ati,2,3 - Solar scale set by \%kk 

For the remaining examples, the bilinear lepton number couplings take the values, 

^ = % = % = L47MeV > ( 4 - 16 ) 

which reproduce correctly the atmospheric mass squared difference, as shown in Fig. 4. A 
single, further lepton number violating coupling \ikk{= ~^kik) is then varied. The branching 
ratio for lepton flavour violating decays when this coupling correctly generates the observed 
solar flavour oscillation are calculated numerically using a Fortran code. 4 We find that there 
are scenarios of this form which correctly reproduce all neutrino data and give rise to branching 
ratios for I — > £'7 which are, or will be, observable in experimental studies. 

When /ii,2,3, Ajfc/c # then diagrams shown in Fig. 5 are generated. It can be seen that 
the amount of mixing on the fermion line inside the loop, again, corresponds directly to 
the amount of mixing between interaction state neutrinos and gauginos/higgsinos. It is this 
4 The same code was used in Ref. [13] to calculate the neutrino masses. 
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mixing which determines the mass of the neutrino produced at tree level and is determined 
by the values given to the bilinear lepton number violating parameters, ^. The left hand 
vertex is determined by the A coupling from the superpotential, as defined in Eq. (1.1). This 
term in the superpotential generates both couplings in the second diagram of Fig. 5. In fact, 
if Xikk ¥= 0, that is, any A coupling with the final two indices the same, a single A coupling 
will generate this diagram. The branching ratio is approximately given by the following 
expression, 



In the following sections, we shall consider in turn all A couplings with symmetric final indices. 

4.2.1 /xi >2 ,3 and A 2 n 

The first example considered is A211. As A 2 n is varied, Fig. 6 shows the resulting solar mass 
squared difference, given by the dashed line, and the \x — > ey branching ratio, given by the 
full line. The light grey strip shows the current experimental value for the solar mass squared 
difference and the dark grey area is the area presently excluded by \x — > ey searches. 

The value of A211 required to generate the correct value for the solar mass difference must 
be comparatively large to compensate for the smallness of the mass of the electron induced 
in the loop. The /j, — > cy diagram generated has a large mass in the loop, and there is no 
suppression in the slepton part of the graph due to intergenerational, or left-right, mixing. 
As such, it can be seen that this scenario, although correctly explaining neutrino masses, is 
ruled out as it predicts a \x — > ey branching ratio which would have been observed already. 
Furthermore, it is shown in Ref. [33,34], that a A211 coupling of this magnitude would violate 
charged current universality. 

4.2.2 /xi >2 ,3 and X 122 

In the top right plot of Fig. 6, it is shown that the value of A122 required to correctly generate 
the solar mass squared difference is smaller; a muon is now produced in the loop contributing 
to the neutrino mass. The lower value of X122, in turn, makes the \i — > ey branching ratio 
produced lower than the previous example, however it is still at the edge of the region ruled out 
by experiment with 90% confidence level. For scenarios with slightly heavier scalar masses 
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Figure 6: \i — > cy and Am^ against A211 or A122 - the dotted line (and right hand axis) 
show the magnitude of Am^ ol with the light grey band indicating the current 3a allowed range; 
the full line (and left hand axis) show the branching ratio of \x — > with the dark grey area 
showing the values currently excluded at 90% confidence level. For the upper two plots, the 
R-parity conserving parameters are set by the SPSla benchmark points, for which the mass 
of the charged scalar that consists mostly of interaction state Jar is approximately 143 GeV. 
For the lower two plots, the mass of the scalars was increased such that the mass of jiR was 
raised to 145 GeV (bottom left) and 265 GeV (bottom right). 
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than the SPSla benchmark point this would not be ruled out. In the two bottom plots 
of Fig. 6 the mass of the scalars has been increased, such that the mass of the charged 
scalar which is mostly fiR is 143 GeV (bottom left) and 265 GeV (bottom right) compared 
to approximately 145 GeV which is produced by the SPSla benchmark values of R-parity 
conserving parameters. We note how sensitive the resulting branching ratio is to the mass of 
the scalar in the loop. As shown in Eq. (4.17), the Br(Z — > I'j) ~ l/m~. Because this scenario 
is at the edge of current limits, and because of this sensitivity to the value of scalar masses, 
it is a particularly interesting scenario which can be studied in future experiments. 

4.2.3 /xi, 2 ,3 and A 3 n,i 3 3 

In the first plot of Fig. 7, we note the effect on the branching ratio of r — > by varying A 311 . 
The values for A311 which correctly reproduce the neutrino data are not ruled out by current 
rare decay searches. Not only is the experimental bound less stringent, but the branching 
ratio is suppressed by a factor of (ra^/m^)(Br(/x — >■ e^z> e )/Br(r — >■ eu T i> e )) ~ 1600 compared 
to the that of Sect. 4.2.1 and the fact that the A coupling itself is smaller. The predicted 
branching ratio generated by the Ai 33 coupling that correctly generates the solar mass squared 
difference, is even smaller due to the lower value of the coupling. 

4.2.4 £ti, 2 ,3 and A 32 2,233 

The right hand plots of Fig. 7 demonstrate that the values of A322 or A233 which reproduce 
the neutrino results are not excluded and well below current experimental sensitivity. We 
note that the values for A322 which produce the neutrino mass are smaller than for A211 due 
to mass of the n produced in the loop and A233 smaller still. 

4.3 Atmospheric scale set by ^1,2,3 - Solar scale set by \' ikk 

When fii, \' ikk # the diagrams shown in Fig. 8 are generated. In this case, the mixing on the 
external leg is driven by the //j term, again being determined such that the atmospheric mass 
difference is produced correctly at tree level. In a similar fashion to the previous section, the 
A' coupling on the left hand vertex is varied and the resulting solar mass squared difference 
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Figure 7: r — > ej (full line, left axis) and Am^ (dotted line, right axis) against Asn^ or 
A133; t — * jjf) (full line, left axis) and Ara^, (dotted line, right axis) against A322 or A233 
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Figure 8: Feynman diagrams generated when \' ikk generates neutrino masses radiatively and 
^1,2,3 induce mixing on the external leg. 
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Figure 9: r — > fij and AmJ oi against A' 122 . 



considered. The branching ratio is approximately given by 
3 
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where .M-, and yW 2 7o are diagonal entries in the squark mass matrices and A4 2 r Tr , is the 

a L art u a .Lit 

off-diagonal term which determines the mixing between the scalar partners of the left and 
right handed quarks. We note that the mixing between e# and charginos is much smaller 
than the mixing of v and neutralinos, and as such there is a suppression relative to the A- 
driven diagrams. Furthermore, at the SPSla benchmark point, the squarks are heavier than 
the charged sleptons; the branching ratios are highly sensitive to the scalar mass and this 
further suppresses A' contributions in comparison with A diagrams. The result being that the 
A' vertices produce a negligible effect in this scenario. 

With only the ^ # 0, setting the atmospheric mass scale, the diagram shown in Fig. 3 
gives a contribution to the r — > fij branching ratio of the order 1.38 x 10~ 20 (Fig. 4) compared 
to which we can ignore the contribution from quark loops, as shown in Fig. 9 which takes A' 122 
as an example. 
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Figure 10: Feynman diagrams generated when the bilinear terms in the superpotential and 
supersymmetry breaking terms, fa and Bj respectively, are non-zero 




16 




Figure 12: Feynman diagrams which contribute to k — > Iff, in the case where trilinear lepton 
number violating couplings are dominant. 



4.4 Atmospheric scale set by ^1,2,3 - Solar scale set by B{ 

The diagrams shown in Fig. 10 are generated when fa, Bj ^ 0. In the first diagram of Fig. 10, 
the mixing on the internal fermion is set by Lij and the mixing on the scalar line is set by 
Bi. Again, we note that the mass of the particle inside the loop for the rare decay diagram is 
of the same order of magnitude as of the particle in the radiative correction to the neutrino 
mass. The contribution to the branching ratio is approximately given by, 

2 

T(k - IjUiVj) . (4.19) 

As such, we see that these diagrams are not ruled out by current experimental bounds, as 
shown in Fig. 11, and are not within reach of upcoming studies. 



T(h - l 3l ) 



I \ 12 „2 

I A oji I e 



(4tt) 2 Glmf c. 




4.5 Atmospheric scale set by A - Solar scale set by A' 

In the following sections, both the atmospheric and solar mass scales are set by radiative 
corrections. Again, we can find combinations of parameters which correctly describe the 
neutrino sector and also give rise to experimentally attainable branching ratios for / — > l lr f. 
We first consider the case in which Ai 33 and A233 are varied over the following range, 

A 133 = -V2A 233 = 5 x 10" 5 - 8 x 10" 5 . (4.20) 

The ratio ensures the correct mixing between neutrino interaction states is reproduced and the 
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magnitude sets in atmospheric masses squared difference. In addition to this, the contribution 
of the first diagram in Fig. 12 to the branching ratio of I — > £'7 is approximately given by, 



T(h - l n ) 



1 \ 12 1 \ 12 r 

\\kk\ \ A jkk\ 



1 1 

24^1 



2 



r(Z< — Zji/iPj) . (4.21) 



(4vr) 2 Gf. 

The results are given in the upper left panel on Fig. 13. The dashed line (and right 
hand axis) show the atmospheric mass squared difference and the light grey band shows the 
values for which Ai33 ; 233 generate an atmospheric mass difference in agreement with current 
experimental observations. The full line (and left hand axis) show the corresponding branching 
ratio for fj, — > e^y. The resulting branching ratio is well below current or future experimental 
sensitivity. 

To generate the solar mass squared difference, X\kk,2kk,3kk are varied in the following hier- 
archy, ensuring the resulting mixing matrix takes the form observed by experiment, 

\l A 2kk A 3fefc ( A r>r)\ 

\kk - - ■ (4-22) 

The second diagram in Fig. 12, produces a contribution to the U — > ^7 branching ratio of 
approximately, 



2 1 v 12 



IV |2 IV 

I ikk I I jkk\ 



1 1 

3 ml 



2 



T(h IjUiVj) . (4.23) 



(4tt)2 G\ 

The results for the three possible cases, that is k = 1,2,3, are shown in the remaining 
panels of Fig. 13. We note that while k = 2, 3 are well below current or planned experimental 
sensitivity, the scenario in which A' m 211 311 generate the solar mass squared difference would 
be discernible by upcoming experimental studies, although we note that bounds from fj, — e 
conversion in nuclei already strongly constrain this set of parameters [35]. When the solar 
mass squared difference is generated by A' 133 233 333 , bottom right panel of Fig. 12, the branching 
ratio is of the same order as that given by A133 233, setting the solar scale, and contributes 
with opposite sign. Resulting in the negative gradient shown. 

4.6 Atmospheric scale set by A' - Solar scale set by 

Again, both the atmospheric and solar mass scales are set by radiative corrections. First, we 
consider the scenario in which the atmospheric mass squared difference is set by A' 111211311 . 
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Figure 13: ji — > e'j anc? Am^ f atm against - In the top left panel Ai.33^33 are varied and the 
resulting branching ratio for /i — > (full line, left hand axis) and Am^ (m (dashed line, right 
hand axis, grey band showing current 3a band from experiment) are noted. In the remaining 
three panels, X' lkk 2kk 3kk are varied and Br(\i — > &~)) and Am^ are shown. 
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Figure 14: \x — > e7 and Ara^ (m against A' m - A' m 211 311 are varied and the resulting branching 
ratio for \i — > e7 f/tt/Z line, Ze/t /iand <ms, dar& gre?/ area showing values excluded by exper- 
imental searches) and Am^ (m (dashed line, right hand axis, light grey band showing current 
3a band from experiment) are noted. 

The parameters are given by, 

V,„ = ^ = ^ ■ (4.24, 

The resulting atmospheric mass squared difference and resulting branching ratio for fj, — > 
e7 are shown in Fig. 14, from which it can be seen that the parameter space which brings 
about the correct atmospheric mass difference is already ruled out by the rare decay searches. 

In a similar fashion, we consider the scenario in which the atmospheric mass squared 
difference is set by A' 133)233i333 , as follows, 

W _ ^233 _ K33 (A 9 r^ 

a 133 - 7^ - 75 • ( 4 - 25 ) 

The results are given in the upper left panel of Fig. 15. In this case, we note that the values 
°f ^'133,233,333 which give the correct value for the neutrino mass, generate negligible rates for 
H -»> erf. 

The solar mass squared difference must now be generated. It can be set either by Ai 33j2 33 
or by a different set of A' couplings. First we examine \' lkk 2kk 3kk , which are varied as follows, 
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The upper right and bottom left panels show the results for these parameters. Again, we note 
that although the scenario in which the solar scale is set by A' 122i 222,322 (bottom left panel) is 
not within experimental sensitivity, A' m 211 311 (top right panel) is close to the current bounds 
and would be seen by searches planned for the near future. 

The final possibility is that the solar mass squared difference be set by A133 233. The 
parameters are varied in the following hierarchy, 

A 133 = -V2A 233 • (4.27) 

The results, given in the bottom right panel of Fig. 15, show that this scenario is well below 
both current and future experimental bounds. 



5 Benchmarks 

Benchmark scenarios for studying R-parity violating models have been suggested [36], for 
which the mass spectrum, nature of the lightest supersymmetric particle and decays have 
been studied. The benchmarks presented in [36] are selected as they produce interesting 
signatures at future colliders and are constrained by measurements of (g — 2) M , the b — > 57 
decay and mass bounds from direct particle searches. 

In addition to these benchmarks, we suggest some other interesting scenarios. Our moti- 
vation being that neutrino data is known and we demand that the model correctly reproduces 
these results. As there are numerous combinations of lepton number violating parameters 
which satisfy this requirement, we consider scenarios for which the upcoming I — > £'7 searches 
will constrain the Lagrangian parameters. 

The following combinations of lepton number violating parameters are considered, where 
all R-parity conserving parameters are set at the SPSla benchmark point, 

• Benchmark Scenario 1 - 

^ = -^L = J^= = 1.47MeV , A 211 = 7.4 x 1(T 4 
v 2 V3 



• Benchmark Scenario 2 



^2 ^3 , \ or,.. 1A -2 



fIl = J-± = ^-L = 1.47MeV , A311 = 3.7 x 10 
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Figure 15: \x — > cy and Am^ atm against A^ - In the top left panel A' 133 233 333 are varied and 
the resulting branching ratio for fj, —> cy (full line, left hand axis) and Am^ m (dashed line, 
right hand axis, grey band showing current 3a band from experiment) are noted. In the top 
right and bottom left panels, A' lllj211311 and A' 122i 222,322 are varied, respectively, and the \i — > cy 
and Amj ol are shown. In the bottom right hand panel A133 233 are varied and the /i — > ey and 
Am^ oZ are given. 
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Benchmark Scenario 3 - 



A 133 = -V2A 233 = 6.5 x 10 - 5 , A' m = ^ = = 1-05 x 

• Benchmark Scenario 4 - 

A' = ^233 = ^333 = 3. 25 x 1Q -5 y = = _^3U = ^ x 1Q -2 

In the first two benchmark scenarios the observed flavour oscillations of atmospheric neu- 
trinos are driven by the bilinear lepton number violating term in the superpotential giving 
rise to a mass difference at tree level. For this to occur, the /i, parameters are of the order 1 
MeV. In Benchmark Scenario 1, the solar mass squared difference is then generated by setting 



ei 



A211 = 7.4 x 10 . Merely for comparison, we note that this is approximately A211 ~ 25 y, 
where y e is the Yukawa coupling associated with a given particle, in this case being the elec- 
tron. This will give rise to branching ratios for /x — > e^y which can be probed by upcoming 
experimental studies. In Benchmark 2 the solar mass squared difference is determined by 
A311 = 3.7 x 10~ 2 ~ 1200 y e . This combination of parameters will generate a branching ratio 
for r — > e7 which may be probed by future studies of r decays. 

In the Benchmark Scenarios 3 and 4 both the atmospheric and solar mass squared differ- 
ences are set by radiative corrections. The bilinear lepton number violating terms are set to 
zero, and the neutrinos are all massless at tree-level. In Benchmark Scenario 3, |Ai 33)233 | ~ 3y e 
set the atmospheric mass difference and |A' m 2 ii 311I ~ 40^ sets the solar mass squared 
difference. In Benchmark Scenario 4, |A' 133 233 333 | ~ 0.1 yd sets the atmospheric scale and 
|A' m 211311I ~ 40 yd gives the solar scale. In both Benchmarks 3 and 4, \x — > cy would give 
branching ratios which will be observed by future experimental studies. 



6 Conclusions 

That lepton flavour violating decays of charged particles have not been observed is worthy of 
note. The suppressed branching ratios arise automatically in the Standard Model, but this 
is not the case in Supersymmetric extensions where it already puts strong bounds on certain 
parameters. In this paper, we have examined the effects of lepton number violating couplings 
on these branching ratios. Combinations of parameters which describe the neutrino sector 
were chosen to be examined and their subsequent effect on the rare decays considered. 
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We first investigated the case in which R-parity conserving parameters are set to the 
SPSla benchmark point and the bilinear lepton number violating term in the superpotential, 
Hi generates the atmospheric mass difference. We showed that the values of [ii which correctly 
describe the atmospheric mass difference and mixing angles are not ruled out by the current 
bounds onr^ /ry, ejov/j,^- e^y. As such, the bounds from lepton decays are less stringent 
than the bounds from neutrino data. We then considered the case in which a further lepton 
number violating parameter correctly reproduces the solar mass difference and considered the 
combined effect on the I — > £'7 decays. We note that in this scenario these decays can impose 
constraints on one of the trilinear lepton number violating parameters in the superpotential, 
A. We considered all the examples in which the A coupling has symmetric final indices, which 
generate the solar neutrino mass with just one non-zero coupling. A 2 n and A122 are excluded 
by experimental searches for \i — > &y; X311, X133, A332, A233 are not. We note, however, that 
the branching ratios are sensitive to the masses of the scalar particles in the loop. As such, in 
scenarios where the scalar masses are heavier than those in SPSla, the branching ratios can 
be greatly suppressed. 

In this scenario, the limits on / — > Z'7 do not place useful constraints on A', or the bilinear 
lepton number violating terms in the supersymmetry breaking part of the Lagrangian, 
generally, the current constraints from the neutrino sector are stronger. 

Second, we considered the scenario in which trilinear lepton number violating couplings 
are dominant. We set all bilinear lepton number violating couplings to zero, and again set Im- 
parity conserving parameters to the SPSla benchmark point. In this scenario, both neutrino 
mass scales are determined by radiative corrections and in order to generate the correct mixing 
matrix in the lepton sector, more than one lepton number violating coupling must be non- 
zero for each mass scale. Because of this, diagrams which contribute to / — > £'7 are generated. 

We note that limits already exist when A' m 2U 311 are used to generate mass differences in 

(') 

this scenario. In general however, for A^,^ the constraints from the neutrino masses are 
stringent. 
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A Appendix 



In this section, we present the Feynman Rules for the J£~MSSM. The mixing matrices, Z, 
diagonalise the mass matrices of the model and determine the amount of interaction eigenstate 
in each mass eigenstate. The full mass matrices and definitions of the mixing matrices are 
presented in Ref. [13]. 

A.l Neutral Scalar - Charged Fermion - Charged Fermion interac- 
tions 



H° q 



K , 



%e 1 ry ry* ry K 1 7 V 7* 



sw a/2 



e 1 

j=Z AlqZ*_ lp Z +2r 
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e 1 
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A. 2 Charged Scalar - Neutral Fermion - Charged Fermion interac- 
tions 
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Charged Fermion - Quark interactions 



a Rk 



+i^'aij Z uiq Z -(2+a)p Z d R jk 



a Lj 



u Lk 



-i — Z tiq Z +lq Z d L ij + i (Xu)i k Z t(Z + k)q Z d L ij Z 



+i (Yu)^ Zl iq Zl Rjk Z +2p 



+i^aij Z d(3+j)q Z -(2+a)p Z u L ik ~ i — Z diq Z -lp Z UL 
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